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ABSTRACT 

Background quasars are potentially sensitive probes of galactic outflows provided that one 
can determine the origin of the absorbing material since both gaseous disks and strong bipolar 
outflows can contribute to the absorption cross-section. Using a dozen quasars passing near 
spectroscopically identified galaxies at z ^ 0.1, we find that the azimuthal orientation of the 
quasar sight-lines with strong Mg II absorption (with W^'^^^^ > 0.3 A) is bi-modal: about half 
the Mg II sight-lines are aligned with the major axis and the other half are within a = 30° of 
the minor axis, suggesting that bipolar outflows can contribute to the Mg II cross-section. This 
bi-modality is also present in the instantaneous star- formation rates (SFRs) of the hosts. For 
the sight-lines aligned along the minor axis, a simple bi-conical wind model is indeed able 
to reproduce the observed Mg II kinematics and the Mg II dependence with impact parameter 
b, (W^'^^^^oc b~^). Using our wind model, we can directly extract key wind properties such 
as the de-projected outflow speed Vout of the cool material traced by Mg ii and the outflow 
rates Mout- The outflow speeds V^ut are found to be 150-300 kms~^, i.e. of the order of the 
circular velocity, and smaller than the escape velocity by a factor of ~ 2. The outflow rates 
Mout are typically two to three times the instantaneous SFRs. Our results demonstrate how 
background quasars can be used to measure wind properties with high precision. 

Key words: galaxies: evolution, galaxies: formation, galaxies: haloes, galaxies: intergalactic 
medium, galaxies: kinematics and dynamics, quasars: absorption lines 



1 INTRODUCTION 

In spite of our understanding of the growth of dark-matter struc- 
tures from the initial density fluctuations (e.g. [White & Rees|1978[ 
|Mo & W hite "2002 ), the halo mass function over-predicts the ob- 
served number density of galaxies both at the low- and high-mass 



ends of the mass function (e.g. 


Croton et al.|2006 Jvan den Bosch 


etal.|2007| 


Conroy & Wechsler 


2009||Behroozi et al.'2010j,Moster 


et al.||2010 


IGuo et al.||2010l |] 


Firmani & Avila-Reese||2010| and 



references therein). This major discrepancy requires a (or several) 
mechanism to somehow suppress galaxy formation. 

Super-novae (SN) driven winds are often invoked because they 
could suppress star formation in low-mass galaxies (L < L*), 
(e.g.|Dekel & Silk|1986 , Oppenheimer & Dave 2006l|Oppenheimerl 
|et al.|2010| and transport large amounts of energy and gas out of 



young galaxies and enrich the inter-galactic medium (IGM). This 
scenario is supported by the fact that galactic winds are ubiqui- 
tous in all types of star-forming galaxies: in local starburst galaxies 
(e.g. Lehnert & Heckman 1996, Heckman et al. 2000, Str ickland] 
& Stevens 2000, Strickland et al. 2004, Martin 1998, M artin et al.| 
,2002^. Schwartz et al.,2006j , in extreme starbursts, such as the Ultra 
Luminous Infra-Red galaxies (ULIRGs) (e.g. |Martin|2005| [Rupke| 



let al.|2 005; 'Martin'20061 [Martin & B ouche 2009|), and in normal 



star-forming galaxies both at intermediate (Sato et al.|2009||Weiner 
et al.|2 009, Rubin et al. 2010a' b]) and high-red shifts ( [Pettini eTaL 



2002c||Shapley et al.,2003,,Genzel et al.|2011| ). 



Numerical simulations must often invoke strong galactic out- 
flows in order to reproduce the luminosity function and the enrich- 
ment of the IGM (e.g. Oppenheimer & Dave|2006| |Oppenheimer| 



let al.|2010|[Schaye et al.|2010[|Wiersma et al.|201 These simula- 
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tions must, however, postulate scaling relations for the wind speeds 
and outflow rates, etc., in order to reproduce observational con- 
straints. For instance, Dave and collaborators assume that the out- 
flow rate, Mout, is proportional to the SFR (Mout = rj SFR), where 
the loading factor 77 is a function of halo mass. For momentum- (or 
energy-) driven winds, 77 is proportional to V~-^ (or V~^), respec- 
tively. 

Unfortunately, most wind properties (e.g. the opening an- 
gle, the outflow rates, the loading factors) are poorly constrained. 
The best estimates of Mout made by several groups over the past 
decades (e.g.peckman et al. 1990 2000, Pettini et al. 2002b , Mar] 
|tin et al.||2002y ^Martin^ 2005) using galaxy absorption line spec- 
troscopy are usually uncertain by orders of magnitude. One reason 
for these large uncertainties is that one must estimate the total gas 
column in the wind from the ion column density which requires 
assumptions for the gas metallicity and the ionization factor. An- 
other reason is that traditional spectroscopy (e.g.'Leh nert & Heck] 
man|1996[|Heckman et all2 000, Martin 1998 Rup ke et al.|2005[ 



Martin|2006[ [Schwartz e tal. 2006, |Weiner et al.|2009TlRubin et al.| 



2010a|b|) probes the wind looking 'down-the-barrel', i.e. it provides 



no information on the physical location of the material, as the blue- 
shifted material could be located at 0.1, 1 or 10 kpc from the host. 
In addition, the wind geometry is unknown and as a result the wind 
solid angle is often assumed to be Qw = 47r. Lastly, the absorp- 
tion trough could be filled with difficult-to-remove emission (e.g. 
[Prochaska et al.||2011| ). Each of these factors make estimates of 
outflow rates very uncertain. 

The radial extent of the wind can be addressed directly with 
background galaxies, as demonstrated by |Steidel et al.| p010| ) at 
z ^ 2, and by Bordoloi et al.|{|2011|) at z ~ 1. However, apart 
from exceptional cases (e.g. Rubin et aL]|2010a| , one must stack 
the spectra of dozens or hundreds of background galaxies in order 
to gain sufficient signal-to-noise. This stacking inevitably leads to 
averages in the geometries involved |Steidel et al.|2010| ). But when 
sufficiently large samples are available, the azimuthal dependence 
can be revealed, as demonstrated by [Bordoloi et al.| ( [201 1| ) who 
showed that the rest-frame Mg II equivalent width of background 
galaxies is strongest along the minor axis. 

One can also use background quasars to probe the radial ex- 
tent of the wind. Background quasars have several advantages to 
the other techniques. For instance, they allow us to probe gaseous 
material of any distant star-forming galaxy irrespective of its lu- 
minos ity (as in^Stocke et al. 2004[ [Bo wen et al. 2005, ,Tripp et al.[ 
|2005[ [Zych et al.|2007| ). Compared to background galaxies, back- 
ground quasars are better probes because they require no stacking, 
i.e. the geometry of the absorbing flow is preserved and no loss of 
information occurs in azimuthal averages. Compared to galaxy ab- 
sorption line spectroscopy, they allow us to probe the material at a 
known distance from the original source. 

The low ionization Mg II doublet (AA2796,2803) seen in back- 
ground quasars (QSOs) is ideal for probing galactic winds as it can 
be observed from z ^ 0. 1 to z ^ 2.2 in the optical and has been as- 
sociated mostly with star-forming galaxies since the work o f|Berg-[ 
'([1988 ), Berge ron & Boissel ( [T99T] ), [Steidel & S argent] ( [1992| ), 



Steidel et al.| ( |1994 ). Unfortunately, the physical origin of strong 



absorbers is still debated. Indeed, Mg II absorbers could p robe the 
K) material entrained in galactic winds (e.g. Nulsen 



cool (T ^ 10 ^ 

jet al.|1998|[Schaye|2001||Martin|2006[pielouche & Bowe n|2010|), 
the outskirts of gaseous disks (e .g. [Prochaska & Wolfe 1997 Stei-[ 



del et al. 2002, Kacprzak et al. 2010 ), the halos of galaxies (e.g. 
Bahcall & Spitzer 1969 , Mo & Miralda-Escude^l996||Lanzetta &[ 



falling material (e.g. [Tinker & Chen|2008| [Kacprzak erari[20T0l 
Stewart et al.||2011|) or a combination of these mechanisms (e.g. 
Charlton & ChurchilT] |1998| ). While mounting evidence points to 
galactic winds for strong Mg II s ystems (jBond et al. 200 1[ [Bouchej 



jet al.|2006| [Menard et al.|2011[ [Nestor et al.|201lt , a direct link 
between low-ionization metal lines and galactic winds has yet to be 
established. 

The debate on the origin of strong Mg II absorbers arises be- 
cause it is difficult to build large samples of individual quasar- 
galaxy pairs. Indeed, at low redshifts (< 0.1), the frequency of such 
pairs is low (e.g. Bowen et al.|1995[[Stocke et aL 2004 , Bo wen et aL] 
12005 ), and at high-redshifts it is time consuming to identify the 
galaxies associated with QSO absorption lines. Even though sig- 
nificant samples of quasar-galaxy pairs are available (e.g.'Chur chillj 
e^al. 1996, Chen et al. 2001 , Noterdaeme et al. 2010, Chen eTal] 
[2010b||Rao et al.|2011||Lovegrove &^Simcoe|2011i ), there are only 
19 z ~ 0.5-1 galaxy-quasar pairs where the host galaxy kinemat- 



|2002|[Emson et al.|2003 


Chen et al.'2005 , Kacprzak et al.'2010). 


excluding the 14 pairs of 


louche et al. ( 2007 ) whose analysis is in 



Bowen[[1992[ [Sternberg et al.|[2002[|Maller & Bullock|2004| ), in 



progress. 

In this paper, we use the unique sample of about a dozen z ~ 
galaxy-Mg II pairs from [ Barton & Cookej ( [2009| ) and [Kacprzak] 
jet al.[ f2011a) in order to investigate the relative orientations of the 
quasar lines-of-sight with respect to the host galaxy orientation. 
Section J2| summarizes the properties of the sample. We show in 
Section|3]that the sample is made of two classes of Mg II absorbers. 
In Section|4] we discuss the physical properties of galactic outflows 
for the sub-sample of pairs related to outflows. Finally, we discuss 
the implications of our results in Section[5] 

Throughout, we use a '737' cosmology, with h = 0.7, = 
0.3, and Qa = 0.7. 



2 SAMPLE 

A recent increase in blue sensitivity of the Keck Low Resolution 
Imaging Spectrograph (LRIS) opened a new redshift window and 
allows the detection of the Mgll doublet down to 2; ^ 0.1 as 
demonstrated by [Barton & Cooke] p009| ). In combination with the 
spectroscopic completeness of the Sloan Digital Sky Survey (SDSS 
r^rk et al.|2000l [Schn eider et al. 2010 ) at that same redshift, one 
has the possibility to study large unbiased samples of nearby Mg II- 
galaxy pairs. [Barton & Cooke[ ( [2009] ) constructed such a sample 
designed to probe for the presence (or absence) of Mg II absorp- 
tion in a well-understood, volume-limited spectroscopic survey of 
galaxies at 2; 0.11 with luminosity Mr < -20.5 (-- L*). Out of 
the 20 sight-lines passing within ISh^l kpc of z ^ 0.1 luminous 
galaxies, six exhibit strong (W^,^^^^^ > 0.3 A) Mg II absorption at 
the same redshift as the galaxy. Kacprza k et al.[ ( [201 la] ) extended 
this sample to 13 such galaxy-Mg II pairs using the same observa- 
tional strategy. 

This is the largest sample of quasar-galaxy pairs at 2; 0, 
which has the advantage that follow-up observations are either 
available in the SDSS database or easy to obtain using 4-m or 8- 
m class telescopes. This strategy led [Kacprzak et al.[ ( [2011a| l to 
measure the rotation curves of the hosts using the Apache Point 
Observatory (APO) and present a detailed comparison between the 
host galaxy kinematics and the absorption kinematics. Hence, the 
galaxies in this sample have reliable systemic redshifts, a key as- 
pect for our study. 

The galaxy SFRs are computed using the Ha luminosities 
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3 RESULTS 



3.1 Two classes of absorbers 



Azimuth Angle (a) 



Figure 1. Schematic diagram of an inclined disk, showing the relative az- 
imuth angle a measured with respect to the galaxy major- axis. 
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Figure 2. The bimodal distribution of azimuth angle \ol\ for our sample 
made of 10 galaxy-quasar pairs where a could be determined. About half 
the sample has \a\ less than 20 deg., and the other half has |q;| > 60 deg. 
The distribution suggests that extended gaseous disks and galactic winds 
contribute significantly to the Mg ii cross-section in roughly equal propor- 
tions. The lack of galaxy-quasar pairs with \a\ between 20 and 60 deg. is 
significant at the > 3.0(7 level (see text). 



[ig a^Sal peter] ( |1955| ) initial mass function (IMF) 
ddening]^ Due to the small angular size of the 



measured from SDSS spectra using the formalism of Kewley & Do- 
|pita|p002 ) assuming 
and no intrinsic reddening 
SDSS fiber apertures, the SFRs were scaled by the ratio of the r- 
band galaxy total counts to those contained within the SDSS fiber 
The sample of 13 galaxy-quasar pairs of [Kacprzak et al.| 
( [2011a| l is made of 11 quasars and 13 galaxies, two of which are 
associated with the same quasar line-of- sight. In this study, we kept 
only the galaxy with the smallest impact parameter, and were left 
with 1 1 unique galaxy-quasar p airs. Table[T]lists the observed prop- 
erties of the sample, taken from |Kacprzak et aL] ( |2011a| ). 



As mentioned in the introduction, the physical mechanism (e.g. 
galactic winds, inf ailing gas, extended disks) that can produce 
Mgll absorption has been debated for decades. As discussed in 
[Kacprzak et al.||2011b ), if there are preferred kinematics and spa- 
tial distributions of Mg II absorbing gas relative to the host galaxies, 
then the absorption strengths would follow a predictable behavior 
as a function of galaxy orientation and/or relative orientation to the 
QSO sight-line, for instance. 

Here, we investigate the distribution of the azimuth angle a 
of the quasar line-of-sight with respect to the galaxy major axis, 
as illustrated in Figure [l] By examining the azimuth angle distri- 
bution, one could conclude that galactic winds dominate the Mg II 
cross-section if the quasars are preferentially aligned with the mi- 
nor axis {a ^ 90), since this is the only mechanism that can pro- 
duce cool material systematically along the minor axis. Alterna- 
tively, one could conclude that the extended parts of gaseous disks, 
or inf ailing material (according to Stewart et al. 201 T) would domi- 
nate the Mg II cross-section if the quasars are preferentially aligned 
with the galaxy orbital plane (e.g. |Charlton & Churchill|1996| . Fur- 
thermore, if there is no preferred a, then this would point towards 
the halo model where the gas clouds traced by Mg II are uniformly 
distributed (e.g. [Lanzetta & Bowen|1992[|Steidel et al.llQOilpo] 
|& Miralda-Escude|1996[|Tinker & Chen|2008| ). 

Looking at the relative distribution of QSO lines-of-sight has 
previously not been possible because of the low number of unbi- 
ased pairs whose kinematic axis is known. Indeed, to our knowl- 
edge, at z there exists only one such sample, the sample of 1 1 
galaxy-quasar pairs of [Kacprzak et al.|(|2011a|), originally from the 



volume-limited survey of^Barton & Cooke[([2609)^ 



We remeasured the azimuth angles because we noticed some 
inconsistencies in the a's reported by [Kacprzak et al.[ ( [2011a| ). To 
do so, we measured the inclinations (i), position angles (PA) of the 
major axis and azimuth angles {a) of the quasar using two methods: 
a visual inspection of the images and a parametric Sersic ( 1963 ) fit 
to the SDSS postage stamp images. The inclination (via its axis ra- 
tio b/a), the Sersic index, n, and the galaxy PA were fitted with cus- 
tom routines, where the Sersic profile is convolved with the image 
Point Spread Function (PSF). We note that our inclination measure- 
ments obtained from visual inspection ('Manual'), from our Sersic 
fits ('Fit') and from a full bulge-disk decomposition by Kacprzak] 
[et al.[ ( [2011a| ) agree well with each other Table |2] lists the inclina- 
tion and azimuth angle measurements. 

We found that the inclinations can be reliably determined, as 
our fitted values are within 10% to those derived by Kacprzak et aL[ 
( [2011aj . On the other hand, we found that our azimuth angles dif- 
fer significantly, and we attribute the difference to a mistake in the 
image orientation. We note that for one galaxy (J161940G1), its 
PA is poorly determined as it is observed almost perfectly face-on 
{b/a ^ 1.0), i.e. its major-axis PA and its azimuth angle is unde- 
fined. Hence, we are left with a sample of 10 galaxy-quasar pairs 
with reliable azimuth angles. 

Figure [2] shows the distribution of \a\. The distribution is 



This is equivalent to assuming a Chabrier (2003^ IMF and At/ = 
0.8 mag of extinction, typical for star-forming galaxies. 
^ At intermediate redshifts, there is the sample of 19 2; ~ 0.5-1 galaxy- 
quasar pairs from ^Steidel et al.[2002[[Enison et al. [2003 [[Chen et al.[2005l 
[Kacprzak et al.[2010| and the 14 pairs discovered by [Bouche et al.[^2007| . 
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strongly bi-modal, with a subset made of four galaxies with small 
|a|'s and another made of six with \a\ ^ 90. In other words, 
the distribution shows that all the quasar-galaxy pairs are either 
nearly aligned with the major or with the minor axis. The lack of 
quasar-galaxy pairs with a between 20 and 60 deg is not consis- 
tent with small number statistics. Indeed, the probability of having 
no quasar-galaxy with \a\ between ^20 and ^60 deg is 0.2% us- 
ing 10^ simulated samples (with A^pairs = 10) drawn from a uni- 
form distribution [7(0, 90). Thus, the central gap in the distribution 
shown in Fig.|2]is significant at the >3.0-cr level. 

In short, we found a strong azimuthal dependence in the pres- 
ence of Mg II, in good agreement with the results of B ordoloi et al.| 
( |2011| ) obtained at z ^ 1 in the Cosmological evolution Survey 
(COSMOS). In stacked spectra of background galaxies, [Bordoloi | 
|et al.| ( [20TT |) also found a strong azimuthal dependence of the total 
(2796A,2803A) rest-frame equivalent width Wr for pairs within 
40 kpc of inclined disks. Similarly, in stacked spectra of thousands 
of local star-forming galaxies from SDSS/DR7, |Chen et al.|p010c| ) 
showed that the blue- shifted Na I D absorption is stronger within 
60° of the minor axis. 

Contrary to stacked spectra, the background quasars (Fig. [2| 
reveal that there are two distinct populations of Mg II absorbers, in 
roughly equal proportions. One population of absorbers could be 
associated with extended gaseous disks (or other processes aligned 
with the disk) and the other could be associated with galactic winds. 
At this stage, the wind scenario is only based on a plausibility argu- 
ment, namely the expectation that all the other physical processes 
(accretion, extended gaseous disks, halo gas) are expected to pro- 
duce absorption that are not aligned with the minor axis. In the 
next section, we will investigate the wind scenario further and at- 
tempt to match the observed Mg II kinematics using a simple wind 
model. This exercise will show that the classification based on |a| 
is strongly supported by our kinematic model. Table|2]lists the clas- 
sifications. 

Having identified two classes of absorbers, we searched for 
differences in the galaxy properties between the two sub-classes. 
In Figure [3] we show the galaxy postage stamp images and the 
quasar location relative to the galaxy major and minor axes. The 
red and blue squares show the relative QSO positions for the sub- 
sample with \a\ > 60 and \a\ < 20, respectively. We find no 
significant differences in the galaxy colors between the two sub- 
samples. The most significant difference between the two is found 
in their SFR distribution. Figure|4]shows that galaxies with high \a\ 
(classified as 'windy') have higher SFRs, while pairs with low |a| 's 
('disky') have lower SFRs. Note the kinematic modeling presented 
in section [T2J] will show that two galaxies are misclassified based 
solely on the azimuth angle criterion. 



3.2 The Mg II kinematics of the Wind sub-sample: Wind 
modeling 

In light of the results presented in section |3.1| we postulate that 
the Mg II absorption is produced by material entrained in galac- 
tic winds for the sub- sample made of galaxy-quasar pairs aligned 
along the galaxy minor axis (with |a| ^ 90). To test this hy- 
pothesis, we construct a simple bi-conical wind model aimed at 
reproducing the observed Mgll kinematics. Our simple galactic 
wind model is made of 10^ 'clouds' distributed in a cone within 
an opening-angle Oma^ (and corresponding solid angle = 



10' 



10 



i 10° 
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Figure 4. SFR as a function of |q;| shows that galaxies associated with bi- 
conical galactic winds (high a) have higher SFRs than those with low a's. 
SFR are computed from the nebular emission lines in the SDSS spectra and 
scaled to take into account the SDSS fiber size. The kinematic modeling 
in section [32l] show that two galaxies (J092300G1 and J114518G1) are 
misclassified based solely on the azimuth angle a criterion. 



maxy [3 The discrete clouds populate the cone from a minimum 



radius Rmin = 5 kpc to a maximum radius i?max ^100 kpc, cov- 
ering the range of impact parameters. We assume that the clouds 
are entrained in the bi-conical wind and are moving at a constant 
velocity (Kut), which is the only model parameter fitted against 
the data. 

All the parameters related to the geometry of the wind can 
be determined from the data. The wind opening angle is ^max ^ 
30 deg according to the distribution of azimuth angles (Fig. |2]) 
since no QSO-galaxy pairs are found beyond ±30 deg from the 
minor axis. The corresponding solid angle is thus of order unity 
with — 0.86. Similarly, the relative geometric orientation of 
the wind with respect to the quasar line-of-sight is also given by 
the data. For all galaxy-quasar pairs, the galaxy inclination (i) and 
the relative orientation (a) of the sight-line are set by the data. The 
only degree of freedom left is to choose whether the cone inter- 
cepted by the quasar sight-line is pointing either towards or away 
the observer. We adopt the convention that x, ^ are the coordinates 
in the plane of the sky, with x along the galaxy major axis and 
consequently, the 2;-axis along the quasar line-of-sight. 



3.2.1 Notes on individual cases 

Figure [5|a) shows an example of the wind model for the galaxy 
J081420G1 (towards the quasar SDSSJ081420.19+383408.3) 
whose inclination is i ^ 35°, and azimuth angle is |a| ^ 80°. 
The top left panel shows the cone view face-on and the top right 
panel shows a side view of the cone. The solid-blue oval represents 
the inclined disk and the black circles represent the conical outflow. 
The bottom left shows the average 2;-velocities of the clouds as a 
function of position. The QSO location is represented as the filled 



^ After some experimentation, we found that a hollow cone with ^min ~ 
10° performed somewhat better in reproducing the shape of the Mgii ab- 
sorption profile. Note the kinematic results (absorption centroid and width) 
are completely insensitive to ^min • 
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(b) 

Figure 3. Position of the quasar line-of-sight distribution relative to the galaxy major axis (x-axis) and minor-axis (y-axis) for the galaxy-quasar pairs with 
high (a, circle symbols) and low (b, square symbols) azimuth angles a. The size of the symbols is proportional to the rest-frame W^"^^^^. The postage 
stamp images show the color SDSS image of the associated galaxy. Blue point source objects with squares indicate the QSO location when visible within the 
field-of-view. 
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circle. The bottom right panel shows the line-of- sight velocity dis- 
tribution of the clouds at the location of the quasar. The distribution 
is convolved with an instrumental resolution of ~ 150 km s~^, cor- 
responding to the LRIS data of |Kacprzak"eraLl ( [20TTat . The LRIS 
spectra is shown in Figure [5jb). The pink (grey shaded) and blue 
(grey hatched) areas show the range of velocities accounted by our 
wind model and by the halo-disk model of Kacprzak et al.|(2011a| ), 
respectively. 

In this particular case, we choose the model where the cone 
is pointing away from the observer since the observed Mg II ab- 
sorption is redshifted with respect to the galaxy systemic veloc- 
ity (0 kms~^). The wind speed Vout is tuned to match the ob- 
served velocity range. For this particular case, we find that Vout 
200 it 25 km produces a good match to the data. 

Figure [6] shows the wind model for J091119G1 towards the 
quasar SDSSJ091119.16+031152.9. This sight-line has one with 
the highest impact parameter (6 = 71 kpc). While the equiv- 
alent width is relatively small (Vl^^^^'^^^ 0.8 A), the absorp- 
tion is spread over a very wide range of velocities (from —300 to 
300 km s~^). In addition, the galaxy is seen almost perfectly edge- 
on (i ^ 80°). As a result, our constant Vout model requires a large 
wind speed of ^ 500 kms~^. We also note that the Mgll profile 
shows little or no absorption around Vsys, which indicates either the 
absorbing material traces the edges of the cone or that this sight-line 
caught the wind as it stalls, in which case our assumption of pure 
radial wind velocities break down. 

Figure [T] shows the wind model for J102819G1 towards the 
quasar SDSSJ102847.00+391800.4. The main difference with the 
previous example is that this galaxy is less inclined with i ^ 55°. 
As a result, the wind produces absorption only redward of the 
galaxy systemic velocity. Because of the large impact parameter 
(b — 89 kpc) of the quasar, we extended our wind model 
to i^max = 140 kpc and the inferred radial wind speed is 
Vout — 300 kms~^. Note that the observed velocity profile is 
rather well reproduced. 

Figure |8] shows the wind model for J111850G1 to- 
wards SDSSJl 14518.47+451601.4. This galaxy has SFR= 
3.75 Mq yr~^, is less inclined with i ^ 30°, and the quasar impact 
parameter is 6 :^ 25 kpc. As a result, this sight-line likely intercepts 
the wind and the disk at the same time, resulting in significant ab- 
sorption at Vsys- As a result, our model, which does not include a 
disk component, cannot reproduce the entire Mg II profile and the 
wind speed is less constrained. In fact, the disk model of Kacprzak 
|et al.| ( |2011a| ) can reproduce the kinematics not accounted by our 
wind model (blue hatched). 

Figure [9] shows the wind model for J225036G1 towards 
SDSSJ225036.72+000759.4. This galaxy has a moderate SFR of 
1.36 Mq yr~^, is highly inclined with i ^ 70° and the quasar im- 
pact parameter is b ^ 50 kpc. The radial wind speed Vout is 
inferred to be ^ 200 kms~^. Given this geometric configuration 
with the inclination approaching 90°, the sight-line probes a wide 
range of velocities, a feature that is very consistent with the ob- 
served Mg II kinematics. Note also that the profile suggests that the 
cone is hollow, i.e. that the Mg II absorption arise on the edges on 
the cone, which is the reason why we used ^min = 15° in this case. 

We next discuss two pairs (J092300G1, J114518G1) whose 
classification based on a was misleading once the global aspects 
of the geometry are taken into account. In particular, the sight-line 
towards SDSSJ092300.67+075 108.2 has the lowest impact param- 
eter (6=12 kpc). Figure [To|reveals that the low impact parameter 
combined with the high galaxy inclination implies that the sight- 
line can intercept the other parts of the galaxy in spite of having a 



high azimuth angle \a\ = 82°. Hence, this galaxy-quasar pair is 
classified as 'ambiguous'. Furthermore, the galaxy has a low SFR 
of 0.02 Mq yr~^, red colors and an early type morphology (Sersic 
index of n ^ 5). 

The other ambiguous case is J114518G1 towards 
SDSSJl 14518.47+451601.4. This sight-line has a low a of 
^ 15°. However, the wind modeling (Figure [TT} reveals that 
the sight-line can very well intercept the wind given the low 
galaxy inclination (i ^ 30°) if the cone extends a bit beyond our 
canonical 30 deg to 40 deg. The wind model can account for the 
Mg II kinematics with an outflow speed of Vout — 125 km s~^. In 
addition, the galaxy SFR is high with SFR= 2.59M0 yr"\ blue 
colors and an exponential profile, which are additional reasons to 
consider this pair as a 'wind' pair. 

3.2.2 Summary of the Wind Modeling 

In summary, we performed simple kinematic models for the quasar- 
galaxy pairs whose azimuth angles |a| are close to 90°. Figures[5}j9] 
show that a simple geometric model can account surprisingly well 
for the observed Mg II kinematics of four out of five QSO-galaxy 
pairs classified as 'windy' and for two other pairs initially classified 
as 'disky' based on their a measurement (see table|2]). Conversely, 
the wind model is unable to account for the kinematics of the sight- 
lines with a ^ 0° since they do not intercept the conical flow. 

In addition, the disk-halo model oflKacprzak et al.| ( |2011a| ), 
shown as the hatched areas in Figures[5}|9] has serious difficulties in 
reproducing the Mgll kinematics for the 'wind' sub-sample since 
the projected disk velocities are approximately km along the 
minor axis. However, in the case of J111850G1, our wind model 
cannot reproduce the entire Mgll profile but the disk model of 
[Kacprzak et al.| f201 la) can account for the extra component. 

Quantitatively, our analysis reveals wind speeds Vout that are 
typically 100-300 kms~^ (listed in table [2]), i.e. are on the order 
of the maximum circular speed Vma^ determined from the rotation 
curves by [Kacprzak et al.| ( |2011a| ). Based on these results, we will 
investigate the wind properties in more detail in the next section §|4| 
We discuss the properties of the 'disk' sub-sample in Appendix [A| 

3.3 Radial Dependence 



Using the sub-sample of z ^ 0.1 QSO-galaxy pairs that we asso- 
ciate with galactic outflows, we investigate the radial dependence of 
the Mgll absorption. Figure 12 shows the observed VFr^^^^^ (cyan 
circles) as a function of impact parameter b. Because the W^^'^^^ 
can be skewed due to the various inclination effects, we corrected 



the WS" 



(red circle) to a common inclination ofi = 90°. The 



corrections are calculated by comparing the total number of clouds 
intercepted at the observed i to that number for i = 90°. The most 
significant correction is for J081420G1 (see Fig. [5]) because that 
sight-line intercepts a small fraction of the cone compared to the 
edge-on situation. 

For comparison, we show in Figure [12] the radial dependence 
of the cool halo gas from the z ~ 0.5 sample of QSO-galaxy 
pairs collected and analyzed in [Kacprzak et al.[ { [2007 [ |2011b| l 
(squares). Since, the |a| distribution is also bi-modal for this sam- 
ple ( [Churchill et al.|2012[[Bordoloi & et al.|2012| ), we only show 
those with |a| > 45° and whose uncertainty in a is less than 30° 
(3a) in order to remove ambiguous cases. Figure s hows that the 
2; = 0.1 and 2; = 0.5 data sets share a common ^^^-b relation 
which goes approximately as b~^ . A formal linear fit to the data 
gives 6 = — 1.1 ± 0.5. 
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Figure 5. (a) Top left: Kinematic model of conical wind for J081420G1 viewed on the sky plane where the x-axis corresponds to the major axis and the 
y-axis with the minor axis. The solid oval represents an inclined disk and the black circles represent the conical outflow. Top right: Side view where the z 
coordinate corresponds to the line-of-sight direction, with the observer to the left. Bottom left: Map of the averaged line-of-sight velocities (kms~^) of the 
clouds. Bottom right: Distribution (convolved with the LRIS resolution) of the cloud 2;-velocities at the observed QSO location. The normalization of this 
distribu tion is inversely proportional to the impact parameter b (see Eq.|B5}. (b) Observed Mg ii kinematics from the LRIS spectra taken by |Kacprzak et al.| 
I20T lal, where km s ^ corresponds to the galaxy systemic redshift. The solid line (red) shows a multi-Gaussian fit to the absorption features to help the 
reader identify the Mg ii absorption profiles. The pink (grey shaded) and blue (grey hatched) areas show the range of velocities accounted by our wind model 
and by the halo-disk model of |Kacprzak et al.lpOl la) , respectively. 



Our radial dependence is significantly different that the one 
determined by [Bordoloi et aL| p011| ) around of z ^ 1 inclined 
disks using background galaxies. However, we refrain from any di- 
rect comparison since the background galaxy technique used by 
[Bordoloi et al.| p011| ) does not give a measure of the collective 
absorption of Mgll 'clouds', but gives a measure of the radial de- 
pendence of the covering fraction Cf (r) [^because the background 
sources are naturally extendedj^ 

Given that there is an empirical relation between W^"^^^^ and 
the total column density N^ii ( [Menard & ChelouchelpOOQI ), we 
show on the right y-SLxis of Figure |12| the corresponding column 
density iV. The expected radial dependence of the column density 
N{b) for an optically thin medium, whose density p{r) is geomet- 
rically diluted [p{r) = po (ro/r)^] , is 



N{b) 



I 



dx 



62 +: 



At z ^ 2, the cove ring fraction was found to go as (6) oc b 



Steidel et al, 



j2010 ), close to the theoretical expectation of Cf oc 6 



(1) 



2/3 



by 

for adiabatically expanding clouds moving in a hot medium jMartin &] 
[Bouche|2009| . 

^ We note the interpretation of the radial dependence in [Bord oloi et al.j 
( [201 1) using isothermal profiles from Tinker & Chen ( 2008) is not appro- 
priate here since the model of jTinker & Chen, ^2008 j was tuned to quasar 
absorbers. 



whe re th e integral is performed perpendicularly to the cone 
(Eq. |B5| . Eq. shows that N(b) is expected to go as oc b~^. 
The solid line in Figure [12] shows that this is very good descrip- 
tion of the data. The scatter around the solid is only 0.20 dex. The 
match between the data and the expected radial dependence may 
seem surprising given that Eq.[T]assumes an optically thin medium 
whereas the Wr^^^^^ absorption can be optically thick. This can be 
understood if one realizes that the W^^'^^^ is also proportional to 
the number of sub-compon ents or clouds (e.g. [Bergeron & Boisse] 
[799T|[Churchill et aL|2003[[Chelouche et al.|2()08[ ). 

In summary, Figure [12] shows that, for QSO-galaxy pairs as- 
sociated with galactic outflows, there is a tight correlation between 
Wr^^^^^ and impact parameter b, following the expected b~^ de- 
pendence. The scatter around b~^ is very small, only 0.20 dex for 
this sub- sample. While this anti-correlation has been known for two 
decades (e.g.'Lanzetta & Bowen|1990 ','Steidel'1995l [Bouche et aL[ 
2006^ Chen & Tinker 2008|, the scatter has been previously shown 
to be 0.5 dex dChen et al.|20T0al [Churchill et al.|2012| ). 

[Kacprzak et al.[ ( [201 lb| ) argues that this scatter is a function 
of the host inclination, and [Chen et aL] pOlOa ) argues that this 
scatter is correlated with stellar mass and perhaps also with SFR 
|Chen et aL||2010bl ). In Appendix |a] we discuss the W,^'^^^^-b 
relation for the 'disk' sub-sample and show that the large scatter 
is caused by the mixing of several physical mechanisms in Mg II 
samples, namely extended gaseous disks and galactic winds (see 
also [Churchill et al.|20 12] [Bordoloi & et al.|2012| ). It remains to be 
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Figure 6. a): Kinematic model of conical wind for J091 1 19G1 as in Fig.|5] b): the observed Mg ii kinematics with respect to the systemic velocity as in Fig. [5] 
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Figure 7. a): Kinematic model of conical wind for J102847G1 as in Fig.|5] b): the observed Mg ii kinematics with respect to the systemic velocity as in Fig. [5] 
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Figure 8. a): Kinematic model of conical wind for J111850G1 as in Fig.|5] b): the observed Mgii kinematics with respect to the systemic velocity as in 
Fi g. [5] Given the l ow inclination, this sight-line is likely contaminated by absorption from the disk. Indeed, the blue hatched area shows the disk-halo model 
of |Kacprzak"et^pOTTa) . 



Sky view 



J225036G1 



100 [ 



50 



^ 



-50 



-100 



• QSO 




«Q= 77° 

SFR= 1.36 Mg/yr 

W = 1.08 A 



Side view 




100 -50 50 100 -100 -50 50 100 

Ax (kpc) Az (kpc) 

Projected <VJx,y)> v; distribution along QSO los 

■160 looor 



100 

50- 

Q. 

^ 0- 

<l 

-50- 

= 220 km/s 
= 100 kpc 
-100- = [15°,3O°] 



-100 -50 50 100 
Ax (kpc) 




-400 -200 200 400 
KK,yJ [knn/s] 





Mg 


11 1 , 

II A2796 - 


1 1 1 1 1 1 




1,1, 



-400 -200 200 400 



, 1 , , , 1 


, 1 , , , 1 , 
Mg 


II A2803 " 








1 1 1 1 1 1 




, 1 , 


-400 


-200 


200 
vel (km/s) 


400 



(a) 



(b) 



Figure 9. a): Kinematic model of conical wind for J225036G1G1 as in Fig. [5] b): the observed Mgii kinematics with respect to the systemic velocity as in 
Fig.|5] 
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Kinematic model of conical wind for J092300G1 as in Fig. [5] b): the observed Mgii kinematics with respect to the systemic velocity as in 
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Figure 11. Kinematic model (a) of conical wind for Jl 1451 8G1 compared the observed Mg ii kinematics (b) with respect to the systemic velocity as in Fig. [5] 
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Figure 12. W^"^^^^ as a function of impact parameter b for QSO-galaxy 
pairs classified as 'wind', i.e. with \a\ > 60. The cyan circles represent 
the observed W^"^^^^, and the red circles show the W^"^^^^ corrected for 
inclination effects (see text), for the edge-on case where i = 90°. The 
squares show the z ~ 0.5 QSO-galaxy pairs of ,Kacprzak et al.|{2011a) for 
a similar sub-sample with the azimuth angle a > 45. The expected radial 
dependence (oc b~^) for sight-lines intercepting a cone at 90° is shown by 
the solid line and is a good description of these two data sets. The right y- 
axis shows the corresponding A^jj j using the W^^^^^^-Nyii relation from 
[Menard & Chelo uche ( 2009) (MC09). The radial dependence of A^jji 
supported by the Mg ii doublet ratio shown in the top panel. 



demonstrated whether the stellar mass or SFR dependence applies 
to the 'disk' or 'wind' sub-sample. 



4 EXTRACTING WIND PROPERTIES 

Having established that the Mg II kinematics in QSO-galaxy pairs 
with \a\ ^ 90° are consistent with intercepting entrained material 
in galactic winds, in this section, we focus on the terminal velocity 
(§|4.1|) and the mass outflow rate (§|4.2|) of these galactic winds. 



4.1 Terminal velocity 

The cool gas in galactic winds traced by Mg II may be driven either 
by the kinetic energy of supernova ejecta from the entrainment in 
the hot wind (e.g. [Chevalier & Clegg[T985l [Heckman et al.| 19901 



[Strickland & Stevens |2000|), by mom entum injection from the ra- 
diation pressure fMuiray et al. '2005" '2011 '. [Sharma et al.|2011| ) or 
by cosmic ray pressure ( Everett et al. 2008 and references therein). 
Because these predict different scalings for the wind velocity Vout 
with galaxy mass and SFR, it is important to investigate whether 
Vout varies with other galaxy properties. 

In our sample of galaxies which have SFRs of a few M© yr~^, 
we find that the outflow speeds Vout are typically 100-300 km 
(listed in Table [2| using the modeling presented in section [3.2[ 
These relatively low speeds are of the order of the circular velocity 
Vmax- We note that Vout appears to increase with impact parameter. 
We caution that a larger sample is required in order to put strong 
constraints on Voutip), a function that is directly related to the ac- 
celeration mechanism in the wind (e.g. [Murray et al.|2005[[20lT] ). 




40 60 

b (kpc) 



100 



Figure 13. Radial dependence of the outflow speed Vout with respect to the 
escape velocity Vesc{b). This shows that the cool material proved by Mg ii 
is traveling at speeds < lesc- 



We find no correlation between Vout and SFR. This could ei- 
ther be due (i) to our small range in SFRs or (ii) to our Ha-derived 
SFRs which may not be related to the SFRs that occurred when 
the material was launched given that the travel time to the observed 
impact parameter can be significant. Typically, the travel time is of 
the order of 0.5 Gyr (see Table|2}. 

Our measurements of outflow speeds allow us to address the 
following question: Are these velocities sufficient to expel the gas 
from the galaxy into the IGM or will the gas eventually fall back 
onto the galactic disk? The escape velocity Vesc for an isothermal 
sphere is 



Ve. 



V2[l + ln(i?vir/r) 



(2) 



where Vmax is the maximum circular velocity (a proxy for Fvir) 
and Rvir is the virial radius. Since our galaxies are L* galaxies 
with halo mass around lO^^M© , their virial radius is approximately 
i^vir = Kiax/10if(z) 250 kpc, whcrc H{z) is the Hubble 
constant at redshift z. Using Eq.|2] the escape velocity Vesc is 2.5, 
2.3, 1.8 xvrot at b =10, 50, 100 kpc, respectively. We find that 
most of our galaxies have wind speeds that are about half the escape 
velocity. 

Figure p3] shows the outflow speeds relative to the escape ve- 
locity as a function of impact parameter b for our sample of 6 
sight-lines. We use the individual rotation velocities (Vmax) from 
[Kacprzak et al.| (2011a) and appropriate virial radii in Eq.[2] This 
ratio yout/Ksc is < 1 showing that the cool material probed by 
Mg II is traveling at speeds close to the escape velocity. Interest- 
ingly, the ratio Fout/V^esc is about unity for only two sight-lines, 
which are the two with the largest impact parameter (J102847G1, 
J0911119G1). 



4.2 Outflow rates 



Given that our a distribution combined with the results of Bordoloil 
[et al.[ ( [20TT] ) and [Chen et al.[ ( [2010cj ) clearly demonstrate the pres- 
ence of bi-conical outflows within ^max ^ 30° of the minor axis 
of star-forming galaxies, we can estimate the cold gas mass out- 
flow rate corresponding to such a configuration. The outflow rate 
Mout for a mass-conserved flow moving at a speed Vout, with a to- 
tal solid angle Qw, is Mout(^) = p{r) ^wT^ Vout- In the case of a 
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radial sight-line looking 'down -the-barrel' , the ou tflow rate reduces 

2000 , Martin 2005 ) 



to Mout — ^v^N ro Vout (e.g. |Heckman et al.|: 

where the column is = p(r)dr — p(ro)rQ. Similarly, in the 
case of transverse sight-lines at impact parameter b the outflow rate 
is: 



Mout (6) 



NH{b)bVo ut 



(3) 



as derived in Appendix [B] (Eg . [B8] | . Inserting the numerical values 
for ^max — 30°, we have 



Mout (6) = 0.41 Mo yr 



-1 M 



1.5 30° lOiQcm^ 200kms-^ 25kpc 



(4) 



where /a is the mean atomic weight. 

The only unknown parameter in Eq.|4]is the total gas column 
A^j^j, since the impact parameter b and the wind speeds Vout are 
directly or indirectly constrained by the observations. Fortunately, 
we can estimate the mean H I column using the empirical Mg II- 
Hl relation from Menard & Chelouche ( 2009 ) (see also|Rao e t al.| 
[20061 |Bouche|[2008l |Rao et al.||20TTl ) used already in Figure [12] 
Here, we make the implicit assumption that the relation holds for 
sight-lines associated with galactic winds. This assumption is sup- 
ported by |Bouche| ( 2008 ) who reported that there are two popula- 
tions of absorbers, one following the H i-W^^'^^^ relation (corre- 
sponding to outflows) with a high metallicity (half solar to solar), 
and another with roughly constant H I with a low metallicity (1/30) 
corresponding to typical DLAs/sub-DLAs. A direct determination 
of the gas column H I would require observations with the Cosmic 
Origin Spectrograph (COS). 

The typical outflow rates derived from Eq. [4] are ^1- 
5 Moyr"^ for all of our galaxies. These are the most precise 
outflow rates derived for star-forming galaxies for the cold (T ^ 
10^ K) gas|^ Indeed, the uncertainties in Mout are entirely domi- 
nated by the A/jj j uncertainties. Assuming 0.25 dex uncertainty for 
A^jj J, the relative accuracy for the outflow rate is cr(M) /M ^ 0.5. 
Compared to traditional spectroscopy, where both A^jjj, are 
only known to orders of magnitude (e.g. |Heckman et al.|2000|[Per] 
[tini et al.|2002bl ), our technique to estimate mass outflow rates is a 
leap forward. 

Given that the travel time of the low-ionization gas to the ob- 
served impact parameters is significant, a few 100 Myr (see Ta- 
ble[2}, we refrain from comparing our outflow rate to the instanta- 
neous SFR. In future work, we intend to compare them to the past 
SFR determined from stellar population analysis. Overall, the cold 
gas mass outflow rate Mout seems to be 2 times larger than the 
current SFR. 



5 CONCLUSIONS 

In summary, we find that the azimuthal orientation of quasars with 
Mg II absorbers (W^,^^^^^ > 0.3A) relative to the host galaxy major 
axis is not consistent with being uniform at the 3.1-cr significance 
level. The azimuth angle distribution is bi-modal with about half 
the quasars aligned with the major axis and the other half within 
a — 30° of the minor axis (Fig. [5]). This bi-modal distribution con- 
firms the presence of the azimuthal dependence of low-ionization 
gas around inclined disks as reported by [Bordoloi et al.| ( |2011| ) at 



^ Galactic winds are multi-phase phenomenon with potentially additional 
mass in the other phases not included here. However, most of the mass is 
likely contained in the cold phase while most of the energy in the hot phase. 



z ^ 1 and by [Chen et al."(2010cj at z ~ and is inconsistent with 
the halo modeToflChen & Tinto| ( [20081 ); [Tinker & ChenH2008l ) 
and lChen et aLl ( |2010a| ). 

We associate the sight-lines aligned with the minor axis to 
sight-lines intercepting bi-conical outflows ('wind' sub-sample) 
and those aligned with the major axis with sight-lines intercept- 
ing the outskirts of galaxies ('disk' sub-sample). The dichotomy in 
azimuth angle is also present in the instantaneous SFR (Fig. [4]). 

Using the 'wind' sub-sample, the data show that the outflows 
traced by low-ionization lines such as Mg II have several properties: 

• The bi-modal distribution of the a angle (Fig. [2]) shows that 
the outflows are rather well-collimated, covering a total solid angle 

— 2 accounting for both sides of the cone. 

• The wind speeds Vont inferred from the Mg II absorption kine- 
matics and a bi-conical wind model (Fig. [5}|9]) are of the order of 
the rotation speed. 

• The wind speeds tend to be smaller than (or equal to) the es- 
cape velocity, indicating that the low-ionization gas is not escaping 
the halo. 

• The radial dependence of the Mg II equivalent width follows 
approximately the expected b~^ dependence for pure geometry di- 
lution with a scatter of 0.24 dex (Fig. [12]). 

• The mass outflow rates are about 2 x the current SFR, ranging 
from 1 to 6 M© yr~^ using the wind speed and the empirical rela - 



Menard & Chelouche 



2009| ). 



tion between W^'^^'^^ and A/jj i (^-g- 

Compared to the orders of magnitude uncertainties in the best esti 
mates from galaxy spectroscopy (e.g. [Heckman et al.|2000[[Fettini[ 
let al.|2002a| ), our mass outflow rates are accurate to within ~ 50%, 
where most of the uncertainty lies in the A^jj j factor. 

In Appendix [A| we show that our bi-conical outflows are con- 
sistent with the inclination dependence reported by Kacprzak et aL[ 
(2011 b) if the azimuth angle is taken into account. In particul ar, the 
scatter in the W^'^'^'^^-b relation is reduced (as in Kacprzak et al. 



201 lb I when a correction of the type X oc 1/ cosi is applied to 



the 'disk' sub-sample. On the other hand, this correction increases 
the scatter for the 'wind' sub-sample, as one might have expected 
since this inclination correction is not appropriate in this case. Inter- 
estingly, the W^^'^'^^-b relation appears to be much steeper (oc b~^) 
for the 'disk' sub-sample than for the 'wind' sub-sample (oc b~^) 
(see also |Churchill et al.|2012"] ). 

Our results open a new and promising way to study the phys- 
ical properties of galactic outflows at high-redshifts using quasar 
absorption lines. In the near future with larger samples, we will 
be able to investigate further the properties of galactic outflows. In 
particular, a larger sample will allow us to test whether the loading 
factor Tj = Mout /SFR is a function of circular velocity Vc, as being 
assumed in numerical simulations by [Oppenheimer et alT] ( [2010| ) 
and others. 
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APPENDIX A: THE KINEMATICS OF THE DISK 
SUB-SAMPLE 

In the main body of this paper, we find that a significant fraction of 
the Mg II absorbers with VF^,^2796 ^^^^ 3^ found along 

the minor axis, i.e. are not co-planar with the galaxy host. On the 
other hand,l| Kacprzak et aL] ( |201 lb| ) argued that Mg II absorbers are 
co-planar based on their finding that the disk inclination influences 
the scatter of the W^'^^'^^-h relation for a sample of 2; 0.5 QSO- 
galaxy pairs. Thus, their results appear to be incompatible with our 
conical winds since they both apply to the same range of 

O.5-3A. 

Here, we attempt to reconcile these two results. In particular, 
we return to the 'disk' sub-sample with (|a| < 45°) since an incli- 
nation effect should be present predominantly for this sub-sample. 
For extended gaseous disks, we expect that W^'^^'^^ is related to the 
path length X intercepted by the QSO sight-line. Because there are 
only three pairs in the |Barton & Cooke | ( |2009| ) and |Kacprzak et al.] 
( [201 laj sample that meet the a criteria, we include pairs from the 
z ~ 0.5 |Kacprzak et al.] { [2011b| ) sample using the same criteria 
(|a| < 45°), excluding pairs where the uncertainty in the azimuth 
angle is greater than 30° (3cr). 

Figure Slra) shows the W^'^^^^ as a function of impact pa- 
rameter 6 for the z = 0.1 {z — 0.5) QSO-galaxy pairs shown as 
the cyan circles (squares) respectively. The solid line shows a fidu- 
cial (xh~^ radial dependence. The top panel shows that the residual 
scatter from this relation are ^0.63 dex. 

For this 'disk' sub-sample, we expect that the absorption 
equivalent width W^'^^^^ is related to the path length X intercepted 
by the QSO sight-line. For a simple slab geometry, the path length 
X is expected to be oc 1/ cos i where i is the disk inclination. 
Hence, one would expect that the observed 



(Al) 



Al 



b) 



of inclined disks, as stated in Kacprzak et al. 



where Wr is the equivalent width for face-on disks. Figure 
shows the W^^^'^^-h relation for the equivalent width Wr corrected 
to a face-on inclination using X = Xq/ cos i. We set Xq to 0.5 cor- 
responding to an averaged inclination of 60° given that the average 
inclination is < i >— 57° (see the appendix of |Law et al.|2009| ). 
The residual scatter in the W^^^^^-b relation is significantly re- 
duced to 0.32 dex, as indicated by the top histogram. This factor 
of 2 improvement in the scatter rms shows that cos i pla ys a large 
role for the W^-~" 
( |2011b| ). 

This exercise shows that, in some cases, the absorption is 
co-planar and coupled to the galaxy inclination ( [Kacprzak et al.| 
[201 lb|) . As stated in |Kacprzak et al.| ( [2011 b), the absorbing mate- 
rial could also be tracing the accretion of baryons since according 
to [Stewart et al.| p011| ) such infalling material is predominantly 
coupled to the galaxy angular momentum and might dominate the 
Mg II cross-section. 

Conversely, for the 'wind' sub-sample presented in sec- 
tion |3.3| the same path-length correction should not apply and as 
a consequen ce it should increase the scatter in the W^'^^'^^-b rela- 
tion. Figure 



A2 



shows the -b relation for the 'wind' sub- 

sample uncorrected (left) and corrected for the disk path length 
(right). This figure shows that the scatter increases from 0.24 dex to 
0.4 dex and demonstrates that the X path-length is not applicable to 
this subset of Mg II absorbers. Hence, our bi-polar outflows and the 
planar effects reported by Kacprzak et al. ( 201 lb ) are not inconsis- 
tent with each other and this exercise demonstrates the importance 
of the azimuth angle in interpreting Mg ll-galaxy pairs. 
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Table 1. Summary for galaxy-QSO pairs. 



QSO 


^abs 




^em 




Mr 


»^max ^^^^ t' 


SFRtt 


5 


Ref. 






(?>) 


(4) 


(5) 


(6) 
y^j 


(1) 

y ' J 




(9) 
y^ ) 




SDSSJUU5244.23— 005721.7 


0.13460 


J005244G1 


0.13426 


1.46/1.23 


-21.40 


144 


0.05 


32.4 


Kll 


SDSSJ081420.19+383408.3 


0.09833 


J081420G1 


0.09801 


0.57/0.28 


-20.13 


131 


1.27 


51.1 


Kll 


SDSSJ091119.16+031152.9 


0.09636 


J091119G1 


0.09616 


0.82/0.34 


-20.98 


231 


0.26 


71.2 


Kll 


SDSSJ092300.67+075108.2 


0.10423 


J092300G1 


0.10385 


2.25/1.40 


-21.58 


108 


0.02 


11.9 


Kll 


SDSSJ102847.00+391800.4 


0.11411 


J102847G1 


0.11348 


0.30/0.13 


-20.22 


162 


3.75 


89.8 


Kll 


SDSSJ111850.13-002100.7 


0.13158 


J111850G1 


0.13159 


1.93/1.82 


-20.40 


116 


1.96 


25.1 


Kll 


SDSSJ114518.47 +451601.4 


0.13402 


J114518G1 


0.13389 


1.06/1.07 


-21.21 


162 


2.59 


39.4 


Kll 


SDSSJ114803. 17+565411.5 


0.10433 


J114803G1 


0.10451 


1.59/1.25 


-21.58 


67 


0.15 


29.1 


Kll 


SDSSJ144033.82+044830.9 


0.11307 


J144033G2 


0.11271 


1.18/0.93 


-20.22 


112 




24.9 


Kll 


SDSSJ161940.56+254323.0 


0.12501 


J161940G1 


0.12438 


0.32/0.28 


-21.09 


74 


0.06 


45.7 


Kll 


SDSSJ225036.72+000759.4 


0.14837 


J225036G1 


0.14826 


1.08/1.11 


-21.47 


240 


1.36 


53.9 


Kll 



(1) Quasar name; (2) Mg ii absorption redshift; (3) Galaxy name; (4) Galaxy spectroscopic redshift; (5) Equivalent widths for Mg Ii2796 and 28 03A; (6) 
Absolute magnitude; (7) Observed rotation curve velocity (km s~^); (8) SFR in Mq yr~^ ; (9) Impact parameter in kpc; (10) Reference: Kl 1 is for Kacprzak 
I etal.K2011a) . 



Table 2. VA and \a\ measurements. 



Galaxy 


i 


i 


PA 


\a\ 


Class 




g-r 


V^ut 


Mout 


tw 




Kll 


Fit 


Manual/Fit 


Manual/Fit 




(Mo yr-i) 




(kms-i) 


(Mo yr-i) 


(Myr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 










a > 60 














J081420G1 


40+2 


35+2 


30/18+2 


79/67+2 


Wind 


1.27 


0.57 


175+25 


2.2+1.1 


290 


J091119G1 


82+2 


75+2 


• • 753+1 


65/63+1 


Wind 


0.26 


1.02 


500+100 


6.8+3.4 


140 


J092300G1 


56+2 


41 + 1 


20/20+2 


82/84+2 


Ambig. 


0.02 


1.2 


200(?) 


1.4+0.7 


n.a. 


J102847G1 


54+2 


49+2 


89/90+5 


76/83+5 


Wind 


3.75 


0.57 


300+25 


1.0+0.5 


250 


J111850G1 


30+2 


34+1 


86/85+5 


86/59+5 


Wind 


1.96 


0.8 


175+80 


6.0+3.0 


140 


J225036G1 


70+2 


69+3 


56/65+1 


77/69+1 


Wind 


1.36 


1.1 


225+50 


2.2+1.1 


250 








q; < 20 
















J005244G1 


45+5 


42+2 


40/43+2 


20/22+2 


'Disk' 


0.05 


1.0 


n.a. 


n.a. 


n.a. 


J114518G1 


34+2 


34+1 


44/38+5 


15/21+5 


Ambig. 


2.59 


0.67 


125+25 


n.a. 


n.a. 


J114803G1 


45+3 


39+1 


27/31+2 


8/10+2 


'Disk' 


0.15 


1.04 


n.a. 


n.a. 


n.a. 


J144033G2 


55+5 


45+2 


69/75+4 


7/13+4 


'Disk' 




0.50 


n.a. 


n.a. 


n.a. 


J161940G1 


12+12 


5+20 


7/-64+65 


unconstr. 


n.a. 




1.0 


n.a. 


n.a. 


n.a. 



(1) Galaxy name; (2) Galaxy inclination i (degrees) from Kacprzak et al. ( 201 la I who used a bulge+disk decomposition; (3) Galaxy inclination i (degrees) 
from a one component Sersic fit; (4) Galaxy position angle (PA) (degrees) measured manually or from our 2D fits; (5) Azimuth angle |q;| of the quasar 
location with respect to the galaxy major axis; (6) Classification of the quasar-galaxy pair. 'Wind' refers to sight-lines whose Mg ii kinematics can be 
explained with our model. 'Disk' refers to sight-lines whose Mg ii kinematics are li kely related to some oth er physical process taking place in connection 
with the major axis. (7) Instantaneous SFR in Mo yr~^ derived from Ho; taken from jKacprzak et aL|j201 la^ assuming a Salpeter IMF and no reddening; (8) 
g — r color; (9) Radial outflow speed in km inferred from the Mg ii kinematics; (10) Mass outflow rates in Mo yr~^ derived from Eq.|4] (1 1) Travel 

time in Myr from the galaxy to the observed impact parameter (6/V^ut). 



APPENDIX B: MASS OUTFLOW RATES 

Given that we are using quasar absorption lines to determine mass 
outflow rates Mout for the first time, we show all the steps in de- 
riving the Mout equation used in this paper (Eq.|4]). 

In most general terms, the outflow rate Mout for a fluid mov- 
ing at a velocity V through an area Q is 



Mout(r) 



dA p{r) V • fi, 



(Bl) 



where fi is the normal to the surface. For a cone of opening an- 
gle ^max, the outflow spccd is normal to A in spherical coordi- 
nate, and this reduces to Mout = po ^ut where Q = 
27r(l — cos^max). Bccausc, the gas column density of a radial 
sight-line is = J^^ dr p{r) = po ro for a fluid obeying the 

continuity equation (p(r)r^=const), the outflow rate reduces to the 



trivial equation Mout = N ro Vout^ |Heckman et al. 2000 Martin 
|2005| ). In the case of a conical geometry, with a transverse sight-line 
at impact parameter b, the outflow rate reduces to a similar form 



Mout OC N{b)bVout- 



(B2) 



For a transverse sight-line intercepting the symmetric z-axis 
at 6 = 62 of a cone, the integral in Eq. |B1| is performed on the 
cross-section A of the cone at bz. Using t as the radius on the 
cross-section A, the velocity Vz normal to A is Vz = V • fi = 
Vout z:^^ - Hence, the outflow rate Mout (bz) is dA p{r) Vz, 



Moutibz) 



po ro 



27V tdt 



1 



Po To 27r bz Vo 



out / 

Jo 



dt 



(6i+t2)3/2 



(B3) 
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Figure Al. (a): W^'^^'^^ as a function of impact parameter h for QSO-galaxy pairs classified as 'disk', i.e. with \a\ < 45. (b): Same as (a) with the W^'^'^^^ 
normalized to the disk path length X = Xq / cos i, where i is the galaxy inclination. The top panels show that the scatter is reduced from 0.63 dex to 0.32 dex. 
This shows that such QSO-galaxy pairs near inclined disks with \a\ ~ are intercepting either the extended parts of gaseous disks or, as argued in |Kacprz^ 
|et al.,p011bj , the accretion material that can also be co-planar according to the recent simulations of ,Stewart et al.,p01 1 j . 




Figure A2. (a): W^.^^'^^^ as a function of impact parameter h for QSO-galaxy pairs classified as 'wind', i.e. with \a\ > 45. (b): Same as (a) with the VK/^^'^^^ 
normalized to the disk path length X = Xq / cos i, where i is the galaxy inclination. The top panels show that the scatter is increased from 0.24 dex to 0.40 
dex. This shows that the disk path length X is not appropriate for QSO-galaxy pairs classified as 'wind', as one might have expected. 



where t is bound to a maximum rn 
edge. After integration, we find 



bz tan ^max by the cone 



Mout{bz) = PO To Vout 27t[1 - COS^max] 

max • 



since the opening angle ^max defines the integration range xi — 

6 tan ^max. 

In the most general case, for a transverse sight-line that is off- 
^^^^ seted from the z-axis by by, where xi — y^tan^ Oma^bl — by, and 



The column density N{b) for a transverse sight-line intercept- 
ing the symmetric z-axis dXb — bz is 



Nib) 



Po ro 



dx 



1 



XI 



^^-^ arctan — 
b 



X ^1 



Po ro 



2(9, 



(B5) 



b = \^bl -\- b1, we have the column density N{b) 



N{bz 



PqtI^ ^ V^tan^ O^^Ji - bl 
— - — 2 arctan - 



(B6) 



which reduces to Eq. |B5| when by — 0. 

Combining Eq.|B4|with Eq.|B5| we have that the outflow rate 
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collimated winds 



determined from transverse sight-lines is: 

Mout{h) = ^^^^Kut27r[l -COs6'naax] (B7) 

^ iV(6)6yout|^max (B8) 



This paper has been typeset from a TgX/ MlgK file prepared by the 
author. 
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